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ABSTRACT 


A  theory  is  presented  which  permits  the  study  of  the  effects  of  horizontal  winds 
on  the  dispersion  and  amplitudes  of  acoustic-gravity  waves  in  the  atmosphere. 

It  is  shown  that  the  effective  horizontal  group  velocity  for  a  given  frequency  in 
a  given  normal  mode  depends  on  direction  of  propagation  as  well  as  on  frequency 
and  thai  it  is  not  necessarily  in  the  same  direction  as  the  horizontal  wa^  e  num¬ 
ber  vector.  A  number  of  useful  integral  theorems  are  derived  from  a  variational 
principle  and  one  is  subsequently  applied  to  the  development  of  a  perturbation 
method  for  t'  c  computation  of  wind  effects  on  dispersion.  Application  of  the 
method  to  a  realistic  example  indicates  that  winds  car  .appreciably  alter  the  dis¬ 
persion  of  the  normal  modes  and  that  they  should  be  considered  in  any  quantita¬ 
tive  interpretation  of  experimental  microbarograms. 
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ILLUSTRATION 


Figure  1  Char aeteriatic  frequencies  und  a*2  veraua 

horizontal  wave  number  k  for  a  layer  with  sound 
speed  c^,  wind  velocity  and  scale  height  H,. 

In  (a),  it  is  assumed  that  k  *  ■  (1/3)  kcai  •  ^ 

(b)»  it  is  assumed  that  k*  >(1/3)  ke^ . .  .  . .  15 

2  (a)  ^proximate  variation  of  parameter  -  q^with 
angtdar  frequency  for  the  fundamental  mode.  The 
dots  represent  values  conqmted  using  profiles  of 
kinetio  energy  computed  by  Pfeifer  and  Zarichny. 

(b)  Plot  of  -  q'  versus  period  (2»/a»)  derived 
from  the  griq»  in  (a)  and  the  group  velocity  curves 
computed  by  Pfeffer  and  Zarichny.  Here  q' 
represents  the  derivative  of  q^  widi  respect  to  &>. 

Note  that  both  q,  and  q,  are  negative . 25 

3  The  net  perturbed  group  velocity  as  observed 
on  aline  making  angles  9  «  0*»  90*,  180*,  and 
270*  with  die  eastward  direction  versus  period 

in  seconds  . . . . . . . . . . 26 
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I.  INTRODUCTION 


Previous  attempts  to  explain  the  features  of  the  microbarograms  of  infrasonic 
waves  recorded  following  nuclear  explosions!''^ and  natural  atmospheric  explo¬ 
sions  9  have  for  the  most  part  been  restricted  to  atmospheric  models  with  no 
ambient  winds.  To  what  extent  this  neglect  of  winds  is  justified  is  not  entirely 
clear.  Diamond^O  has  recently  discussed  their  effect  on  the  apparent  sound- 
speed  profile  above  White  Sands  and  has  exhibited  data  which  would  seem  to 
indicate  that  typical  winds  arc  of  sufficient  strength  to  have  an  appreciable  ef¬ 
fect  on  sound  propagation. 

A  noteworthy  beginning  in  the  development  of  a  theory  which  considers  the  ef¬ 
fects  of  winds  was  made  by  Weston  and  vanHulsteyn! !  who  showed  that  the 
linearised  equations  of  hydrodynamics  arc  still  separable  if  the  winds  are  hori¬ 
zontal  and  vary  in  direction  and  magnitude  only  with  altitude.  They  also  indi¬ 
cated  how  one  might  calculate  the  variation  of  the  horizontal  phase  velocity  with 
frequency  for  fixed  direction  of  the  tiorizontal  propagation  vector  k. 

Pridmore-Brown!2  also  derived  the  general  linearized  equations  for  sound 
propagation  in  an  atmosphere  with  arbitrary  sound- speed  profile  and  horizontal- 
wind  profile.  In  some  respects,  his  theory  went  further  than  that  of  Weston  and 
vanHulsteyn,  in  that  it  dealt  with  waves  from  a  point  source  rather  than  with  free 
waves  (whose  wavefronts  are  vertical  planes).  However,  Pridmore-Brown  con¬ 
sidered  only  the  steady  state  case  as  he  was  not  concerned  with  dispersion 
phenomena.  Furthermore,  since  he  was  interested  in  sonic  frequencies  of  the 
order  of  100  cps  (as  opposed  to  infrasonic  frequencies  of  the  order  of  10“^  cps), 
he  was  enabled  to  make  a  number  of  approximations  which  cannot  be  justified 
for  lower  frequencies. 

In  the  present  ptaper,  the  theories  of  Weston  and  vanHulsteyn  and  of  Pridmore- 
Brown  are  extended  to  the  consideration  of  the  propagation  of  infrasonic  waves 
from  an  idealized  point  source  characterized  by  an  arbitrary  time  variation f  (t). 
The  expressions  derived  for  the  pressure  on  the  ground  at  a  large  distance  from 
the  source  represent  an  extension  of  the  method  of  normal  modes  to  the  propa¬ 
gation  of  infrasonic  waves  from  a  point  source  in  the  presence  of  horizontal 
winds.  Our  derivation  of  these  expressions  is  similar  to  mat  of  Pridmore- 
Brown  and  is  therefore  given  as  briefly  as  possible.  One  substantial  departiiri* 
from  Prid.T>ore*Brown'8  method  appears  in  the  method  ot  treating  the  two»!uld 
integration  over  the  components  of  the  horizontal  wave  number.  It  is  ovir  opinion 
that  the  .na'.hematical  approximation  utilized  by  Pridmore-Brown  as  indicated 
by  Rqn.  (Z2)  in  his  paper  is  not  justified.  His  approximation  would  indicate  that 
the  horizontal  wave  vector  point  i  rar'ialty  from  the  source  aod  hence  in  thr  sarm- 
direction  as  the  e-oup  velocity.  The  theory  presented  in  this  {laper  indicates 
that  this  is  not  tht  case  in  general  and  gives  a  method  of  computing  the  angular 
deviation  of  the  horizontal  wave  vector  from  the  direction  -»if  the  horizontal  grvu.p 
velocity. 
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The  implementation  of  the  theory  reels  on  the  solution  of  two  coupled  first 
order  differential  equations,  which  represent  generalisations  of  the  residual 
equations  discussed  by  Eckart.  With  appropriate  boundary  conditions  (whose 
rationale  we  discuss)  these  are  eigenvalue  equations  fc«r  the  magnitude  of  the 
horizontal  propagation  vector  f.  The  eigenvalues  k  will  depend  on  the  direction 
of  k  when  winds  are  included  as  well  as  on  the  frequency  and  the  mode  index  n. 
The  theory  presented  in  sections  IV  and  V  shows  how  the  group  velocity  may  be 
calculated  from  a  knowledge  of  the  partial  derivatives  of  the  eigenvalues  with 
respect  to  frequency  and  angle  of  wave  normal  to  a  given  horizontal  direction. 

In  section  VI  we  introduce  a  variational  technique  which  leads  to  a  number  of 
integral  theorems  relating  the  eigenfunctions  of  the  residual  equations  and  their 
eigenvalues.  In  particular,  the  theorems  give  methods  of  computing  the  partial 
derivatives  of  the  eigenvalues  from  a  knowledge  of  the  eigenfunctions  for  a  single 
choice  of  parameters  and  therefore  eliminates  numerical  differentiation.  An 
integral  expression  for  the  group  velocity  is  then  easily  obtained. 

In  section  VII  we  discuss  the  case  of  propagation  in  an  isothermal  atmosphere 
with  constant  winds.  This  is  a  particularly  simple  case  and  one  which  should 
be  carefully  studied  before  proceeding  on  to  more  complicated  model  atmos* 
pheres.  Our  analysis  shows  that  both  the  surfaces  of  constant  phase  and  con* 
slant  arrival  time  are  circles  whose  centers  move  with  the  wind  velocity  and 
whose  radii  increase  at  the  sound  speed. 

In  section  VIII  we  use  one  of  the  integral  theorems  derived  in  section  VI  to  de¬ 
velop  a  perturbation  method  for  taking  winds  into  consideration.  This  method 
makes  calculations  of  dispersion  effects  of  winds  highly  feasible  an*  requires 
only  that  the  wind  independent  eigenfunctions  be  known.  These,  however,  have 
been  explicitly  or  implicitly  obtained  by  all  previous  writers  who  have  com¬ 
puted  horizontal  phase  and  group  velocities  for  wind  independent  model  atmos¬ 
pheres.  As  an  example,  we  make  use  of  the  computational  results  of  Pfeffer 
and  Zarichny^  io  find  the  effects  of  a  wind  profile  exhibited  by  Diamond^®  on 
the  dispersion  of  the  fundamental  mode. 
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II.  FORMAL  SOLUTION  OF  THE  LINEARIZED 
EQUATIONS  OF  HYDRODYNAMICS 


The  atmosphere  is  assumed  to  be  an  ideal  gas  which  obeys  the  equations  of 
hydrodynamics.  The  ambient  variables  p^,  and  t  are  assumed  to  satisfy 
these  equations  in  the  absence  of  any  external  perturbation.  The  w'nds  are 
assumed  horizontal  and  independent  of  the  horizontal  coordinates  r  and  y  and  of. 
time  t.  The  linearized  equations  of  hydrodynamics  as  derived  by  Pridmore- 
Brown  will  therefore  govern  the  spatial  and  temporal  variations  of  the  first 
order  quantities  p,p,  and  il.  In  a  somewhat  altered  notation,  these  equations 
are 


D,  (fJyii)  +  PjjWdv/dz  =  -  Vp  -  gP?2 

Djp  +  V-(Ppil)  =  4))rf(t)  Sff  -  ?^)  (2.2) 

D^P  -  w(gp„  +  c^dp^/dz)  =  c^D^p  (2.3) 

The  term  on  the  right-hand  side  of  Eqn.  (2.  2)  has  been  included  to  take  into 
account  the  presence  of  a  source  at  the  point?,, .  The  function  f(t)  is  to  be 
chosen  such  that  the  linearized  equations  predict  as  accurately  as  possible  the 
known  properties  of  the  acoustic  field  in  the  near  vicinity  of  the  source.  A 
direct  intepretation  of  f(t)  may  be  obtained  by  integrating  both  sides  of  (2.  2) 
over  the  volume  of  a  small  sphere  of  radius  R.  Upon  applying  the  divergence 
theorem,  one  finds  in  the  limit  of  sufficiently  small  R  that 


n  da  =  4  IT  f(c) 


(^.4) 


where  n  is  the  unit  normal  to  the  surface  >  of  the  sphere  of  radius  R.  The 
integral  on  the  left  may  be  interpreted  as  the  mass  expelled  from  the  interior 
of  the  sphere  per  unit  time.  (This  assumes  that  negligible  mass  has  been  added 
to  the  air  by  the  explosion  itself.) 

A  formal  solution  of  the  linearized  equations  may  be  obtained  by  use  of  Fourier 
transforms.  Using  this  method,  one  finds  the  excess  pressure  at  a  point  r  at 
time  (  to  be  given  by 


P  (f  .w>  j  (...I  d.j 
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p(?.  «» 


(2.*.) 


where  g(<u)  is  the  Fourier  transform  of  f(t)  .  such  that 


(2.6) 


and  where  P(i,  ut)  may  in  turn  be  represented  as  ajtwo>fold  integral  over  the 
components  k^.ky  of  the  horizontal  wave  number  k, 


P  (r,  w)  =  ir~ 


'// 


d^ke 


Z(z,  2^,w,k) 


(2.7) 


The  expression  Z(z, z^,  u,  k)  must  satisfy  a  differential  equation  which  may  be 
derived  in  the  manner  indicated  by  Weston  and  vanHulsteyn  and  by  Pridmore- 
Brown  for  the  quantities  IlCz)  and  x'(z)  in  their  respective  papers.  The  differ¬ 
ential  equation  we  find  by  this  method  is 

I  (^'  ^)  ^^2  -■  A2i|  Z  =  ..(i/ri)S(z-z„)  (2.8) 


where 


Aj2  -  Q2  -  , 

(2.9) 

A2j  -  (k/fl)2  -  c  '2  . 

(2. 10) 

A  -  {l-l/2y)g/c2  -  (2c2rldc2/dz 

(2.11) 

«B  “  (y  -  1)  >  (g/c2)dc2/dz 

(2.  12) 

and  Q-u-k'7.  Both  A  and  (Vaisala-Brunt  frequency)  are  functions  of  z  which 

characterize  the  ambient  atmosphere.  A  somewhat  more  convenient  represen¬ 
tation  of  Eqn.  (2.  8)  is  that  of  two  coupled  first  order  differential  equations,  i.  e. 

dY/dz  -  AY  -  Aji  Z  -  -  (iAi)^t  -  Zo» 

(2.  13a) 

dZ/dz  ♦  AZ  -  AjjY  .  0. 

(2.  13b) 

The  utility  of  taking  the  equations  in  this  form  has  been  demonstrated  by  Eckart  . 
Following  Eckart,  we  shall  refer  to  these  two  coupled  equations  as  the  residual 
equations.  The  second  of  three  may  be  considered  as  the  definition  of  the  auxili¬ 
ary  function  Y  (i,  <».  k)  . 


-4- 


0 


The  delta  function  on  the  right-hand  side  of  £x]ns.  (2.8)  and  (2. 13a)  reflects 
the  presence  of  the  source  at  altitude  and  requires  there  be  a  discontinuity 
in  Y  and  dZ/dz  at  z  =  z^,  . 

Boundary  conditions  on  the  set  of  coupled  equations  (2.  13)  are  chosen  tc  insure 
that  w,  the  vertical  component  of  particle  velocity,  be  zero  at  the  Earth's  sur¬ 
face  and  are  chosen  to  insure  that  the  total  solution  (2.  5)  conforms  to  causality. 
The  former  requirement  leads  to  the  condition  y=0at  z  =  0.  (The^ derivation  of 
this  condition  is  similar  to  Weston  and  vanHulsteyn* s  derivation^*  of  d  x/dz  =  0 
at  z=0 . )  The  causality  requirement  has  bcten  shown  by  the  author^^  to  be  satis¬ 
fied  if  Z(z,  Zq,  k)  for  real  ,  ky  is  an  analytic  function  of  a>  for  all  complex 
M  in  the  upper  half  plane  and  vanishes  as  <•>  approaches  i«> . 

To  relate  this  requirement  to  one  governing  the  behavior  of  Z  and  Yat  large  z, 
it  is  convenient  to  assume  that  the  properties  of  the  uppermost  region  of  the 
atmosphere  are  such  that  the  z-variation  of  Z  and  Y  above  some  height  z„  may 
be  written  dov^'n  explicitly.  For  this  reason,  we  assume  that  the  atmosphere 
is  isothermal  above  z^  and  therefore  has  a  constant  sound  speed  in  the  upper¬ 
most  region.  The  winds  at  this  height  will  also  be  assumed  to  be  constant  in 
magnitude  and  direction. 

With  this  assumption  as  to  the  nature  of  the  upper  atmosphere,  one  may  show 
that  the  causality  requirement  requires  that  Z  be  of  the  form 

Z  =  D  (2.  14) 

for  z  >  z,g  ,  where  D  is  independent  of  z  and 

«  -  {'m  (“i  - 

is  a  function  of  f.,,  k^  ,  and  ky  .  Here  n^=  (y/2)g/c^  ,  =  (>  - 

where  and  represent  the  ambient  sound  speed  and  wind  velocity 
in  the  isothermal  layer. 

To  insure  causality  and  the  desired  behavior  of  Z  as  a  function  of  complex 
one  must  require  that  a  be  an  analytic  function  of  in  the  upper  half  plane  and 
that  it  approach  i^  as  ^approaches  x, .  This  requirement  may  then  be  used  to 
specify  the  phase  of  «  for  real  '‘'and  l(,  giving 


Ph(a)  -  0 

w  ^  (ii2 

(2.  15a) 

-  ir  '  * 

<  W  '  €^2 

(2.  15b) 

•  ir 

Wq  *  4ii  '  fi.'  J 

(^.  I5c) 

.  0 

j  '  •  (00 

(2.  15d) 
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Ph  (a)  a  n/2 


(2. 15e) 


«u_2  <  «u  <  <a_j 

*  IT  <o  <  ^_2  (2.15f) 

where  the  five  characteristic  frequencies  are  given  by  the  follov/ing  expressions 


‘"12'  ^ 

(2. 16a) 

(2. 16b) 

"0  -  k  • 

(2. 16c) 

with 


S  -  (“in  ^  4 (2.17a) 

T  -  (S2  -  (2.17b) 

The  upper  boundary  condition  on  the  solution  of  Eqns.  (2. 13)  is  therefore  that 
Z  be  of  the  form  (2. 14)  when  i;  >  where  the  appropriate  phase  of  a  must  be 
determined  from  Eqns.  (2.15)  and  (2. 16). 


Following  a  method  used  previously  by  Haskell^  we  may  formally  write  the 
solution  for  z(z,  zo>  “•  of  Eqns.  (2. 13)  in  terms  of  quantities  Z/(z).  Y;  (z>  and 
2:^(2),  Yy(z)  which  satisfy  the  homogeneous  equations  (i.  e. ,  Eqns.  (2. 13)  with 
the  omission  of  the  source  term)-  The  set  Z/,  Yf  is  defined  as  satisfying  the 
boundary  condition  Y/  «  0  at  z  -  0  ,  while  the  set  Z^,  Yy  is  defined  as  satisfying 
the  upper  boundary  condition.  In  terms  of  these  quantities,  the  solution  for 
Z  of  the  inhomogeneous  equations  which  satisfies  both  upper  and  lower  boundary 
conditions  is  given  by 


Z  (z,  z^.  <u,  k) 


-iZ„(z>)  Zj(z<) 

0(Zo)W(*o) 


(2.  18) 


where  z>  and  z<  refer  to  the  greater  or  lesser  of  z  and  z^  ,  and  the  Wronskian 
W  (Zq)  is  defined  to  be 

W  (*„)  -  Y^  (.„)  Z  (.„)  -  Z„  (.,)  Y^(.„)  (2.  1 9) 

One  may  shou*  directly  from  the  Homogeneous  form  of  EIqn)>.  (2.  13)  that  the 
Wronskian  is  independent  of  altitude.  Thus,  we  may  set 

•  (!„)  -  »(0)  -  Y^(0)Z^(0)  (2.20) 


-6. 


where  we  have  made  use  of  the  fact  that  (0)  =  0  . 

A  more  convenient  expression  for  P(?,  w)  may  now  be  obtained  by  inserting  the 
expression  above  for  Z  into  Eqn.  (2.7).  Since  observations  of  infrasonic  waves 
are  usually  made  on  the  ground,  we  take  z  =  0  in  the  resulting  expression.  This 
gives 

In 

Here  the  integr^ion  is  expressed  in  cylindrical  coordinates;  \  representing 
the  angle  which  k  makes  with  the  x-axis.  The  magnitude  of  the  horizontal 
projection  of  is  abbreviated  by  Rjand  the  angle  between  this  projection 
and  the  x-axis  is  denoted  by  S . 


0 


k  dk  e 


ikR-i<os«9-^) , 


(2.21) 


ni.  THE  METHOD  OF  NORMAL  MODES 


While,  in  principle,  the  order  of  integration  in  Eqn.  (2.21)  is  immaterial, 
great  care  should  be  exercised  in  choosing  this  order  if  one  or  both  integra¬ 
tions  are  to  be  performed  approximately.  The  technique  utilized  by  Pridmore- 
Brown  was  to  first  integrate  over  ^  using  the  saddle -point  method  and  then  to 
integrate  over  k  using  the  method  of  residues.  In  our  opinion,  tliis  order  of 
integration  leads  to  incorrect  results  since  the  saddle -point  method  is  inappli¬ 
cable  if  the  saddle-point  is  close  to  a  pole.  This  appears  to  be  the  case  in 
Pridmore-Brown’s  method,  for,  in  utlizing  the  method  of  residues  in  the 
integration  over  it  is  the  behavior  of  the  integrand  near  its  poles  which  is 
of  principle  importance.  This  objection  can  be  overcome  if  one  first  does  an 
integration  over  k  using  the  method  of  residues  and  then  does  tlie  integration 
using  the  saddle -point  method.  This  is  the  program  we  follow  here. 

The  path  for  the  k  integration  in  Eqn.  (2.21)  is  first  deformed  to  one  which  (in 
tlie  case  of  u>  '>  o)  encloses  all  poles  of  the  integrand  in  the  first  quadrant  and 
which  includes  a  contour  going  from  the  origin  along  the  positive  imaginary 
axis  as  well  as  contours  around  branch  lines.  (All  branch  lines  emanating  from 
branch  points  in  the  first  quadrant  are  taken  as  extending  vertically  upwards. ) 
For  large  Rj,  the  predominant  contribution  to  the  total  integral  comes  from  the 
residues  of  those  poles  which  lie  on  the  real  axis.  The  remaining  terms  may 
be  discarded.  The  integrand  thus  obtained  for  the  0^  integration  will  consist  of 
a  sum  of  terms,  each  representing  the  contribution  from  one  of  the  poles  in  the 
integration  over  k.  Each  term  is  integrated  separately,  and  the  saddle -point 
approximation  is  assumed  to  be  applicable  in  each  case.  The  saddle -point  is 
taken  as  that  of 

eip  k„  R-J  cos(^^  -  (3.1) 

where  k^  (u,  0^^  )  is  the  location  of  the  nth  pole  in  the  integrand  prior  to  the 
performance  of  the  k  -integration.  (One  should  note  that  the  saddle -point  will 
generally'  not  be  at  -  o.)  In  general,  the  validity  of  the  saddle -print  method 
may  be  expected  to  increase  with  increasing  Rj.  If  no  real  saddle -point  c  xi'^ts 
for  any  particular  term  we  may  assume  the  contribution  from  that  term  to  be 
negligible  for  large  R-j-  relative  to  any  other  term  with  a  real  saddle -point. 

The  resulting  expression  for  :^(f ,..)  appropriate  for  large  Rj  appears  as  a  sum 
over  normal  modes  in  the  form 


n 


-  H- 


Both  the  i^se  factor  ^  and  the  amplitude  factor  A„  for  each  normal  mode  are 
/unctions  of  to  and  6,  but  not  of  Rj . 


The  method  for  obtaining  tiie  quantities  and  /3„  follov/s  from  our  preceding 
remarks  concerning  the  procedure  for  derivii^  £qn.  (3.?.).  For  (o>0,  let 
be  a  real  positive  root  of 

- .0  ''*•  ' 


where  n  =  1,  2,  3, .  • .  is  an  index  distinguishing  the  roots.  The  labelling  is 
chosen  such  that  is  piecewise  continuous  in  ^  aikd  ut.  Then,  let  {&>,  0  ) 
be  the  saddle -point  of  the  expression  (3. 1),  or,  equivalently,  a  root  of 


I (w ,  cos -  6)1  »  0 


(3.4) 


For  weak  winds  (which  is  the  case  of  physical  interest)  it  should  be  required 
that  -  6j  is  leas  than  n/2  in  the  event  that  one  need  distinguish  among 
multimple  roots  of  Eqn,  (3.4) 

In  terms  of  k  =  k„  6^  )  and  (<j,6  ),  one  may  set 


=  k  cos  -  0/ 

A„  -  2(2/o)^/2Q|j^.-i(r/4 

where 

Q’‘  -  ^  in(*o)Yu(0)/7^^^)i 

ff  ,  „  -  t  k|^  (i.) ,  6^;  cos  (6|j  -  6)1 

To  compute  and  a„  for  given  uj  and  6,  one  first  determines  the  appropriate 
values  of  and  k  and  uses  these  in  the  above  expressions.  The  differentia¬ 
tion  in  Eqn.  (3.  7)  is  carried  out  at  constant  ,  while  that  in  ( ?.  8)  is  carried 
out  at  constant  6,  (A  more  convenient  expression  for  '2  is  given  by  Eqn.  (6.  14).) 

It  is  not  necessary  to  consider  the  case  <->  o separately,  since  the  fact  that 
both  f  (« )  ar i  r(>*  • '  )  rc&l  implies 

fit.,.,)  -  (-.9) 


(3.5) 

(3.6) 

(3.7) 

(3.8) 


.9. 


and  therefore  implies 


A„(«.d)  -  An(-<u,d)*  (3.10a) 

.  -/5„(-«)  (3. 1  Ob) 

(It  may  be  assumed  that  is  real,  since  we  are  limiting  ourselves  to  undamp¬ 
ed  modes. ) 


The  svirfaces  of  constant  phase  for  given  m  and  mode  number  q  are  determined 
by  the  condition  that  in  the  exponent  of  Eqn.  ( 3. 2)  be  constant,  or 

R.  i-  <3.1 

where  K  is  a  constant.  The  normal  to  such  a  surface  at  a  given  value  of  6 
makes  an  angle 

d-t»ii“^(Rf^  dKj/dd) 


with  the  x-axis.  It  is  readily  shown  from  Eqns.  (3.4),  (3.5),  and  (3. 11)  that 
this  angle  is  identical  with  ^  )•  Thus  the  vector  k  is  perpendicular  to  the 
surfaces  of  constamt  phase. 


•  10* 


IV.  GROUP  VELOCITY 


The  expreasion  (2. 5)  for  the  preaeure  on  the  grouxul  aa  a  function  of  time  may 
be  rewritten  uaing  the  approximation  (3. 2)  in  the  form 


(4.1) 


where 


Re 


(4.2) 


repreaenta  the  normalized  contribution  to  the  total  waveform  from  the  nth 
normal  mode.  Although  the  integration  limita  are  written  aa  o  and  it  muat 
be  borne  in  mind  that.  gener<'illy«  &e  aolution  to  Eqn.  (3.  3)  correaponding  to 
the  nth  normal  mode  will  exiat  only  for  a  limited  range  of  <u.  There  may  be 
either  an  uppi^r  or  a  lower  cutoff  frequency  •  or  poasibly  both.  It  ia  conceivable 
that  the  mode  may  alao  have  a  nxunber  of  frequency  gapa.  for  which  the  mode 
doe  a  not  exiat.  To  allow  for  auch  aituationa.  we  adopt  the  convention  that  the 
A„  {cj,  B )  ahould  be  conaidered  aa  zero  whenever  there  ia  no  corresponding  root 
of  Eqn.  (3,  3), 

The  traditional  method  of  evaluating  the  integral  in  Eqn.  (4.  2)  ia  the  method  of 
atationary  phaae^^.  Although  the  method  in  ita  unn^odified  form  haa  limited 
applicability  to  acouatic^gravity  wavoa.  the  modificationa  deviaed  by  Scorer*^ 
and  by  Weston^  may  generally  be  incorporated  to  make  the  method  a  valid 
approximation  for  the  computation  of  the  waveform  at  diatancea  greater  than 
$000  km  from  the  aource.  For  manv  qualitative  aapecte  of  the  interpretation 
of  the  waveforma,  the  unmodified  method  appeara  to  be  aatiafactory.  In  thia 
paper  we  restrict  ouraelvea  to  the  traditional  method.  A  direct  applicetion 


E 


(ij#,  )  ew*  {  fiiuJ  i  (4.  S) 


where 


(;«>*  •  Ai.w,  «»1 

(  i  *  ' 


(4.4) 


II- 


(4.5) 


and  aij  =  61- (<),  j/R-j.)  is  a  root  of 
tJ/3  (<u,  (i)/dbt  ^  t/R-j- 

(The  liubscript  n  is  omitted  for  brevity. )  The  sum  in  (4.  3)  extends  over  all 
such  roots  if  more  than  one  exists.  The  quantity  $  is  the  phase  of  g{t.^  and 
the  parameter  f  is  it  /2  or  0,  depending  on  whether  ^  ^  /du>^  is  positive  or 
negative,  respectively. 

The  concept  of  group  velocity  is  derived  from  the  method  of  stationary  j^se. 
The  time  t  obtained  from  Eqc,  (4.  5)  represents  the  time  relative  to  the  excit* 
ation  of  the  sourc  ;  at  which  a  wave  of  frequenc/  m  in  the  nth  normal  mode 
airrives  at  a  point  described  by  the  coordinates  Rj  ,  0  .  Thus,  we  may  consider 
the  magnitude  of  the  horizontal  group  velocity  as  being  giveii  by 

Vg  (c,  ^  *  [d  (4.6) 

Since  we  are  assuming  that  the  medium  does  not  vary  in  the  x  ory  directions, 
it  must  be  assiuned  that  the  group  velocity  is  in  a  radial  direction  away  from 
the  source  (i.  e.,  in  a  direction  making  an  angle  of  ^with  the  x-axis). 

In  terms  vf  the  parameter  which  describes  the  direction  of  the  horizontal 
wave  vector  C  with  respect  to  the  x-axis^  the  components  vg,  and  vgy  of  the 
horizontal  group  velocity  are  given  by  the  expressions: 


«* 


cos  ^  f  (sin 
O  k^/<i  w ) 


(4.  7a) 


V 


sin  ^  -  k"*  (cov  \  )  k^,  ^  ) 

(crk^  do) 


(4.  7b) 


where,  in  evaluating  the  partial  derivatives.  lt„(a*  )  in  a*  found  from  Eqn. 
(3.  3). 


The  proof  of  Eqns.  (4.  7)  follows  from  Equs.  {3,  4)  and  (3.  5),  which  give 

s*e  ~  ft)  ^ 

^  {tv*  (6^  -  ft)  i  f#  d  ti! 

that  the  magnitwfs  of  is  the  same  regardless  of  whether  tt  or 

is  kept  constant  wlule  differentiating. )  hisertion  of  these  expressions  into 
Eqns.  (4.  7)  gives  a  group  velocity  with  rne  magnitude  (4.  6)  and  with  a  direci.ion 
making  an  angle  of  ^  with  the  x«axis. 


n~ 


It  shovild  be  noted  that  surfaces  of  equal  phase  are  not  necessarily  surfaces 
of  equal  arrival  time.  If  a  K  independent  of  0  can  be  chosen  in  Eqn.  (3. 1 1) 
such  that  the  resulting  value  of  1  is  equal  to  {u)t  0  )t  lor  some  time  t , 

this  would  be  the  case.  However,  this  would  require 

d 

—  =  0 

which  would  in  turn  require  that 

d  \  1  I 

— -  -  -  V  ..  0 

96^  j  da  ^ 

which  is  clearly  not  true  in  general.  Two  circumstances  where  the  above 
would  be  satisfied  are:  (a)  no  winds  {k  independent  of  ^);  and  (b),  k(tj,  ) 

directly  proportional  to  o.  The  latter,  as  we  show  in  section  VII,  occurs  for 
the  case  of  an  isothermal  atmosphere  with  constant  winds. 


V.  EIGENFUNCTIONS  AND  EIGENVALUES 


For  all  practical  purpoaes,  the  only  modes  of  inteiest  which  may  oe  attained 
from  Eqn.  (3.  3)  are  those  for  which  (z,  ot,  k,  ii<  *ero  at  z  •  An  additional 
mode  apparently  exists  where  0  which  apparently  represents  a  disturb¬ 

ance  traveling  in  the  direction  of  the  ambient  wind  at  altitude  with  the  same 
velocity  as  the  wind.  One  may  discard  this  mode  if  he  is  interested  only  in 
modes  which  travel  with  speeds  of  the  order  of  the  sound  speed  at  the  ground. 
(We  assume  the  wind  speeds  are  significantly  less  than  the  rpeed  of  sound. ) 

If  k  =  k„{tu, is  a  root  of  the  equation  Y^(0)  =  0,  then  the  corresponding  pair  of 
functions.  Zy  (z)  and  Y„(2),  may  be  considered  as  an  eigenfunction  pair  of  the 
coupled  differential  equations  (or  residual  equations). 

dY  dz  -  AY  -  A,,  Z  -  0  (5.  la) 

dZ/dz  *  AZ  -  Aj2  Y  0  (5.  lb) 

and  may  be  denoted  by  the  symbols  Z^(z,  (^)  and  Y^  (z.  w,  6>)c).  The  root  0^) 

may  be  considered  as  an  eigenvadue.  The  problem  of  finding  the  roots  cf  Eqn. 
(5.3)  may  therefore  be  considered  as  that  of  finding  the  eigenvalues  of  Eqns. 
(5.1). 


In  accordance  with  our  remarks  in  section  III,  we  consider  only  those  eigen¬ 
functions  which  correspond  to  real  eigenvalues.  Thus,  A  j2  and  Ajj  nraust  both 
be  real.  This  implies  that  any  set  of  solutions  of  Eqns.  (5. 1)  which  satisfy  the 
lower  boundary  condition  of  Y  =  oatz  =  0  must  be  real  functions  of  z,  ipart  from 
a  multiplicative  constant  which  may  be  complex.  Thus,  the  upper  boundary 
conditions  cai'.not  be  satisfied  if  the  phase  of  a  in  Eqn.  (2.  14)  is  0  or  » .  The 
phase  of  a  must  be  »  2  .  This  proves  that  any  real  eigenvalue  k(€.>,  satisfies 
the  condition 

W  I  (k,  ^  Oj  '-i/ly  (k.  ) 

whexe  wj  axxj  u>2  af®  given  by  Eqns.  (2.1b).  The  nature  of  such  a  constraint. 

IS  best  demonstrated  by  plotting  the  functions  and  In  figure  1,  these, 
as  well  as  are  plotted  versus  k  for  fixed  angle  between  k  and  r,,  .  Nun-tcrical 
values  used  are  such  that  in  (a),  k*.  -(l  .tc^k,  and  in  (b),  k- ^  -(1  hc,„k  ,  For 

simplicity,  k  is  plotted  in  units  of  J  and  «.<  is  plotted  in  units  o.  Kn.  where 
it  the  scale  height  >n.  Figure  1  may  be  coniiderco  as  representing  general¬ 
isations  of  the  diagnostic  diagram  for  a  quiescent  isothermal  atmosphere  given 
by  Eckart 
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ftiVt  NtJMBt**  k  ,  uBilt  o*  i/M^ 

alj». 

F.-tuf,  1  CHA^AC  It«!ST!C  r«tOULNv‘.tS  ,  ,  A»lD  VlRSl’S  HORIZONTAL 

W<VL  NliMBFIl  .  FOR  A  layer  WITH  SOUND  SFttO  ,  WIND  VELOCITY 
A  »ii  S(  Al.E  >  CICHT  :  IN  (.0,  IT  IS  ASSUMED  that  •  1  3 

while  in  (S  ,  II  IS  ASSUMED  THAT  .  -  U  >!  - 


VI,  INTEGRAL  THEOREMS 


A  number  oi  theorems  may  le  derived  relating  the  variation  of  the  eigenvalues 
and  eigenfunctions  of  Eqns.  (5. 1)  to  small  variations  in  the  atmospheric  sound 
speed  and  wind  speed  profile.)  and  to  small  variat‘:ions  in  tbe  parameters  u  ana 

These  theorems  all  follow  from  a  general  theorem  which  we  now  prove. 

We  assume  that  Y  and  Z  are  a  set  of  eigenfunctions  of  Eqns.  (5.  1)  for  given 
A(z),  A 12  (z),  and  A22  (z)  .  Let  5A,  5Aj2.  and  5A21  be  small  variations  in  the 
quantities  A,  Aj2,  i'Ad  A2] .  Then  let  Z  4  SL,  Y45Ybe  solutions  of  Eqns.  (5. 1) 
when  A,  A]2,  ai^  A21  are  replaced  by  A  4  5A,  A22  4  5A21  and  A22  4  5A22  .  (It  is  not 
n.icessarily  assumed  that  Z  4  5Z,  Y4  £Yare  a  set  of  eignfunctions,  but  it  is  assumed 
tl!'.at  they  conform  to  the  upper  boundary  condition.)  To  first  order,  the  variations 
S7.,  SY  will  then  satisfy  the  two  coupled  inhomogeneous  equations 

d(SY)/di  -  A(SY)  -  A22  (5Z'  -  (SA)Y  4(5A2i)Z 

d(SZ)/dz  4  A(SZ)- Aj2(5Y)  =  -(AA)Z  4  (SAj2) Y 
It  follows  directly  from  Eqns.  (5. 1)  and  (6. 1)  that 

|y(SZ)  -  Z(5Y)|  =  (6Aj2)Y’  -  (SA2i)z2  -  2  (5A)YZ 


(6.  la) 
(6.  lb) 

(6.2) 


The  desired  theorem  is  now  obtained  by  integrating  both  sides  of  the  above 
equation  with  respect  to  z  from  0  to  giving 


•• 

|z(aY)|j,0  =  j  |(5Aj2)Y2  -  (5A2,)Z2  ..  2(aA)Yz|  dz 


(6^3) 


where  use  has  been  made  of  the  fact  that  both  Y  and  Z  approach  zero  as  z 
proaches  infinity  and  of  the  fact  thatY>.0  at  z  -  0. 

Under  the  most  general  variation  we  will  consider,  c^,  v,  k,  u.  and  go  to 
4  5c^,  etc.  The  corresponding  variations  SA,  SAjj  to  first  order  may 

be  found  from  Eqns.  (2.9*2.  lZ)to  be 

5A  -  -(A  c^)«c^  -  (2c^)"‘  (d  dz)  (Ac^) 

AAj.,  -  -  (|/c^>  (d/dz)  (Ac^)  -  2  Ok.  A* 

-  2  «(k.  V  fc)  Ak  .  2  n  A  c,  .  2  0k»  »in  Atl),  (6.  S) 


.16. 


^^21  ~  ^  (2k^/fl^)k.  S  V  J-  2k^  5k 


-(2k2/n3)  _  (2k3v/n3)  sin(^  -  (?^)  5^ 

where  0^  (z)  ia  the  angle  between  v  and  the  x-axis. 


(6.  t) 


Insertion  of  expressions  (6,  4-6.  6)  into  (6.3)  w^ith  an  additional  integration  by 
parts  to  eliminate  terms  with  (d/dz)  8c^  givej 


Q  ~  **  +  I2  ^  ^  ^^k  ‘^1  **  ^2 


where 


(6.7) 


Ij  =  2 


J  |n(k.v/k)Y2  +  kt^n-^z^l  dz 


(6.8) 


‘2  =  2 


ac 

J  jflY^  (k2/n3)z2|  dz 


(6.9) 


I3  =  2k 


^  V  sin  )  |n  Y*  i  (k^  n3)z2|  dz 


\  "3 


J  |[>  r  c-»  -  n2  ,2]  y2  _  ,n,2^2 

c 

.  j^(2<t  ch  (k  0>*  -  (t  c-*]  Yzj-  (5c2) 

J  |uv*’  .  (k^u ')/•’!  ^'5;  J, 


dz 


(6.  10) 


(6.  11) 


In  obtaining  the  expression  for  \j.  use  has  been  made  of  the  differential  equa¬ 
tions  (5.  1)  to  eliminate  terms  in  d/  dz  and  dY  dz  . 
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Wc  may  now  apply  Eqn.  (6.  7)  to  a  number  of  special  cases; 

1.  Expression  for  dY/dk  at  z  ^  0. 

We  consider  SI  as  being  zero.  Then  (6.  7)  gives 

(^Y/dk)^,Q  -  -  Ij/ZW  (6.  13) 

The  factor  Q  which  appears  in  Eqn.  (3.  7)  is  therefore 

Q  =  -  Z  UJ  Z  (0)/  (Q  (2^)  Ij  1  (6.  14) 

2.  Expressions  for  and  dk/d\. 

The  variation  of  k  with  for  a  given  normal  mode  and  for  fixed  may  be 
obtained  from  (6.  7)  by  letting  Thus 

dk/3<o  =  I2/I1  (6.15) 

In  a  similar  manner,  one  finds 

dk/ae^  =  Ij/Ii  (6. 16) 

3.  Expression  for  group  velocity. 

Insertion  of  Eqns.  (6.15)  and  (6.16)  into  Eqns.  (4.7)  gives  expressions 
for  the  components  of  the  group  velocity,  ^ter  some  manipulation,  one 
may  write  the  resulting  two  equations  as  a  si  igle  vector  equation  in  the 
form 


g 

where 


V  -  (!j/l2)kA  +  <7 


f  -  (2/I2) 


/ 


.  {k^'n^)Z^  i  d: 


(6„  17) 


(6.  18) 


One  should  note  that,  in  the  event  *  is  constant,  '  ^  is  « 

4.  Eff'ict  of  Atmospheric  Perturbations  onk  . 

11  the  si  jenfunctions  /  and  Y  and  corresponding,;  eigenvalue  k,.  for  g.ven 
and  u  a:vd  for  a  given  model  atmosphere,  are  known,  then  the  cff'-ct  on  the 
eigenvalue  k  of  varying  the  atmospheric  sound  and  wind  profiles  is  given 
to  first  order  by 
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Sk  =  (Aj  -  Aj)/!! 


(6.19) 


The  consequences  of  this  equation  are  discussed  in  section  VIII. 
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VII.  PROPAGATION  IN  AN  ISOTHERMAL  ATMOSPHERE 
WITH  CONSTANT  WINDS 


The  simplest  model  atmosphere  including  winds  is  one  with  a  consta. it  temper¬ 
ature  and  a  constant  wind  velocity.  For  such  an  atmosphere,  the  upper  layer 
with  constant  sound  speed  c^,  and  wind  velocity  coincides  with  the  entire  at¬ 
mosphere.  The  solutions  of  Eqns.  (5.  1)  which  satisfy  the  upper  boundary  con¬ 
dition  are  of  the  form  (2,  14).  The  lower  boundary  condition  requires 

i  a  =  -  A  -  -  (I  -  y/2)g'c2  (7.  1) 

It  follows  that  there  is  only  one  eigenvalue;  it  being  given  by 


k  =  u>/  [c„  ^  cos(^  - 

The  corresponding  functions  Zand  Y  may  be  taken  as 


Y  =  0 


The  relationship  (3.4)  between  6  and  becomes 

-  ^vm) 

t*n  ^  *  - — - -  - 


(7.2) 


(7.  3a) 
(7.  3b) 


(7.4) 


or  (omitting  the  subscript  m) 

c  sin  -  0)  -  V  sin  ~  6^)  (7.5) 

The  corresponding  expression  for  fiU>.  '/)  is  readily  found  to  be 

_ _ _  _ _ _  (7.6) 

I  I  -  (v/ c)7  sin^  ) !  *  7  ,  (v.'c)  cos  ((>  -  ) 

and  the  group  velocity  is  therefore 

~  ^  7  .  ,  (rf  ~  0^  )  (7.7) 

which  is  independent  of  frequency.  The  significance  of  the  above  equation  is 
much  clearer  if  it  is  expressed  in  terms  of  cos  -  *\)  and 

iin  (H  ~  fty)  With  a  little  manipulation.  Eqn.  (7.  7)  then  assum'es  thr 
form 
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which  shows,  as  might  well  be  expected,  that  the  propagation  is  equivalent  to 
that  from  a  source  moving  with  velocity  v  in  a  medium  with  sound  speed  c. 

The  surfaces  of  equal  arrival  time  are  circles  whose  centers  are  displaced 
from  the  origin  in  the  direction  of  the  ambient  wind.  The  ratio  of  the  radius 
of  any  circle  to  the  distance  of  its  center  from  the  origin  is  c/v  .  Since  /i(w,  rf) 
is  directly  proportional  to  to,  the  surfaces  of  constant  phase  will  coincide  with 
surfaces  of  equal  arrival  time. 
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VIIL  PERTURBATION  MET!IOD 


An  approximate  method  of  incorporating  winds  into  a  theory  of  acoustic -gravity 
wave  propagation  is  a  perturbation  method  based  on  Eqn.  (6.  19).  The  unper¬ 
turbed  atmosphere  is  taken  as  one  in  which  v=o  .  The  eigenvalues  are  then 
approximately  given  (to  first  order  inv)  by 

k(uj.  '•  (w)  —  qj(w)co«  ^  -  qy  (<i))sin  (®'  M 

where  q^  and  qy  are  components  of  the  vector 


1 


(lu) 


y2  i  vdz 


(8.2) 


The  first  ter  n  k®(ti)  in  Eqn.  (8.  1)  is  the  unperturbed  eigenvalue.  The  second 
and  third  terms  are  the  first  order  correction  tok(u.  8)^)  and  are  derived  from 
Eqn.  (t).  19)  with  5v  replaced  by  v.  The  quantities  Zand  Y  in  Eqn.  (8.2)  are 
the  sercth  order  eigenfunctions  and  are  computed  assuming  there  are  no  winds. 
The  cuantities  k^^iu)  .  q2(<u)  .  and  qyfcu)  will  be  independent  of  but  will  depend 
on  the  mode  index  n  as  well  as  <u. 


The  apparent  uncoupling  between  th;  i\  rnd  o  dependences  in  Eqr.  (8.  1)  makes 
the  resulting  formulas  for  (f)  and  the  surfaces  of  constant  phase  relatively 
sim  pie.  To  first  order  in  k®  and  it®,  one  has 

f  Q 

/f  (o.  ^  -  k®  (ai)  cos  rt  —  q^,  (o)  sin 

The  surfaces  oi  constant  phase  (see  Eqn.  (3.  1 1)  )  to  the  same  order  of  approx- 
i  nat)0'.t  are  given  by  the  equation 


-  :Kq^  ‘  .  |y  -  (Kqy  k°  )*  )|  ’  -  (k*’)^ 

(8.4) 

where  K  is  a  constant  and  g  -»  R-r*^***^  .  y  -  Ry  sm  .  The  surfaces  are  there- 
f  >re  Approximately  circles  whose  centers  are  displaced  from  the  origin  in  the 
dit  action  of  o  and  which  are  characterised  by  the  number  i  ^  k®  which  gives 
the  ratio  of  the  distance  of  each  circle's  center  from  the  origin  to  its  radius. 
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The  group  velocity  may  be  readily  computed  from  Kqns.  (4.  6)  and  (8.  3).  To 
first  ordc-  one  finds 


''g  'g  ^  [*** 

where  is  the  zero  order  group  velocity,  or  (<?  k«  /(V  I  while  q^’  and  q^'  are 

the  derivatives  with  respect  to  t,. of  q^  and  q^  . 

The  surfaces  of  constat. t  arrival  time  are  circles  to  the  same  degree  of  approx¬ 
imation 


1^*  -  (v®)^  q;  tj  2  y  -  (v^)‘  qy  tj  ^ 

The  center  of  the  circle  moves  with  a  velocity 
out  from  the  center  with  a  velocity  v®. 


(H.6) 


{v0)2  r^q/do;  while  the  waves  move 


To  demonstrate  the  utility  of  the  method  outlined  above  ,  we  apply  it  to  the 
computation  of  the  effects  of  winds  in  a  realistic  case.  The  wind  profile  is 
taken  as  measured  by  Diamond  ^  above  White  Sands  and  as  exhibited  in  figure 
i  of  his  paper.  The  cited  figure  ind’cate.s  that  north-south  components  may  be 
neglected.  Taking  x  to  be  in  the  eastward  direction,  we  therefore  have  q^,  o  • 


In  the  computation  of  q^  we  rely  on  the  numerical  computations  of  Pfeffer  and 
Zarichny^.  In  particular,  we  use  the  plots  of  kinetic  energy  versus  altitude 
for  the  52  km  model  a.s  given  by  figure  5a  in  their  paper  for  periods  48,  87, 
225,  and  3i5  seconds.  The  kinetic  energy  tiiey  tabulate  should  be  proportional 
(as  regards  variation  with  z  )  to  the  quantity 


./'*  7* 


1., '  Y  ‘ 


in  the  notation  used  here.  For  low  frequencies  and  for  the  52  km  model  atmos¬ 
phere  used  by  Pfeifer  and  Zarichny,  it  appears  that  one  may  safely  neglect  the 
second  term  m  the  above.  Thus  the  quantity  qo.)  is  approximately 


(8.  ') 


Curvesi  ol  and  versus  prficxl  a»e  i;tvrn  in  tiijure  ^  of  PfeffVr  afid 

/arichuv’s  papee. 


Using  Eqn.  (8,  7),  we  have  computed  q(w)  ^  q,  for  the  four  frequencies  for 
which  KE  is  plotted  by  Pfeffer  and  Zarichny.  A  curve  fitted  to  these  values 
is  given  in  figure  Za.  A  good  fit  to  this  curve  over  the  frequency  range  of 
interest  is 

q,  M  -  -2.7  V  i(f~^ 

where  cj  is  in  radians /sec  and  q,  is  in  meter Thus  the  magnitude  of 
increaseis  slowly  with  increasing  frequency. 

In  figure  /lb  we  plot  the  factor  versus  period.  Curves  obtained  using 

Eqn.  (8.  5)  of  group  velocity  versus  period  are  given  for  various  values  of  & 
in  figure  3.  It  should  be  noted  that  the  princ  ipal  effect  of  winds  in  this  partic-^ 
ular  example  is  to  increase  the  group  velicity  for  downwind  propagation  and 
decrease  it  for  upwind  propagation  by  an  increment  of  roughly  15  meter/sec. 
However,  the  winds  also  strongly  affect  the  dispersion  of  tlie  wave.  Since  the 
group  velocity  curve  for  dow'T  wind  propagation  is  flatter  than  that  for  upwind 
propagation,  the  signal  observed  to  tlie  east  (i^jwind)  will  be  more  dispersed 
than  that  observed  to  the  west  (downwind). 
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IX.  r,r\’r.i  nni\’G  REMARKS 


The  theory  outlined  in  this  paper  gives  a  method  whereby  the  effects  of  v  inds 
may  be  readily  incopcrated  into  the  sio:dy  of  the  propagation  of  acoustic -gravity 
waves.  Furthermore,  the  example  treated  in  the  previous  section  indicates 
that  the  consideration  of  winds  may  t  ■  nec  ‘ssary  in  a  quantitative  interpreta¬ 
tion  of  act  lal  microbarograms. 

The  principal  complication  introduced  by  v  inds  is  that  they  transform  the  at¬ 
mosphere  into  an  anisc*Tcpic  medium.  Tne  magnitude  of  the  horizontal  wave 
number  ve  ctor  (which  acts  as  an  index  of  refraction)  depends  on  the  vector's 
direction  ,;is  well  as  on  frequency.  This  would  appear  to  make  the  computation 
of  the  pha;5c  £uid  group  velocities  more  difficult  by  an  order  of  magnitude.  How¬ 
ever,  the  perturbation  theory  developed  here  (which  takes  advantage  of  the  small 
ratio  of  wind  velocity  to  sound  speed)  requires  only  the  computation  of  two  func¬ 
tions  t,,  and  q^.  i-  I  in  addition  to  the  wind -independent  wave  number  (,,■>  . 

With  this  simplification  the  consideration  of  winds  becomes  feasible. 

The  p>e " tuvbation  method  is  but  one  application  of  the  integral  theorems  develi’pvd 
in  this  pai’cr.  These  show  promise  of  being  useful  in  the  numerical  calculation 
of  phase  and  group  velocities  as  well  as  in  the  development  of  approximate 
methods  ?f  solving  the  residual  equations. 

The  question  now  remains  as  to  whether  or  not  a  model  atmosphere  with  winds 
independent  of  horizontal  coordinates  is  a  satisfactory  model  for  the  actual 
atmosphere.  Certainly,  it  should  be  more  satisfactory  than  a  nodel  atmosphere 
without  u  nds.  However,  a  glance  at  the  flow  patterns  of  the  atmospheric  winds 

1  w 

on  a  global  scale  given  in  the  Handbook  of  Geophysics  indicates  that  some 
modifica  ion  of  the  theory  may  be  required  to  take  into  accc.unt  the  curvature  of 
the  strcijnlines  cf  tl  ambient  winds.  Such  a  modification  should  be  necessary 
for  prop  igation  Ov-er  horizontal  paths  of  7,000  km  or  greater.  Th>e  present 
theory  n  ay  be  readily  extended  to  cner  such  situation'^  by  using  the  mathemat¬ 
ical  techniques  discussed  by  the  author*^  in  a  previous  theory  of  wave  propaga¬ 
tion  in  an  aimost-stratified  medium.  This  extension  wi.l  be  given  in  a  later 
article. 
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